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We study the superfluidity of a pure spin current that is a spin current without mass current. We 
examine two types of pure spin currents, planar and circular, in spin-1 Bose gas. For the planar 
current, it is usually unstable, but can be stabilized by the quadratic Zeeman effect. The circular 
current can be generated with spin-orbit coupling. When the spin-orbit coupling strength is weak, 
we find that the circular pure spin current is the ground state of the system and thus a super-flow. 

We discuss the experimental schemes to realize and detect a pure spin current. 
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I. INTRODUCTION 

Since the experimental realization of Bose-Einstein 
condensation in optical traps, much effort has been de¬ 
voted to the study of spinor superfluid EH3- With the 
degree of freedom of spin, spinor superfluid has much 
richer phases than scalar superfluid as it has both su¬ 
perfluid order and spin textures. However, most spinor 
superfluids studied till now carry both spin current and 
mass current. It will be interesting to see whether these 
two currents can decouple and further whether a pure 
spin current which has no mass current can flow friction- 
lessly. Our motivation also originates from condensed 
matter physics, in which the concept of spin supercon¬ 
ductor [8, j9[, formed by the Cooper-like pairs of electrons 
and holes and carrying spin super-current, is proposed. 
It is interesting to have these ideas realized in the field 
of cold atoms. 

In this work we focus on the unpolarized spin-1 Bose- 
Einstein condensate (BEC), where the pure spin current 
can be generated by applying a small magnetic gradient. 
It was found in Ref. [IJJ] that a planar pure spin current 
in such a system is always unstable as the m = 1,-1 com¬ 
ponents can collide into the m = 0 component destroying 
the spin current. We find that the pure spin current can 
be stabilized with the quadratic Zeeman effect and be¬ 
come a super-flow. Furthermore, we propose a scheme 
to create pure spin current at the ground state, thus free 
from the issue of instability. Our scheme utilizes the spin- 
orbit coupling. Specifically, we study a spin-1 BEC with 
Rashba spin-orbit coupling confined in a two-dimensional 
harmonic trap, and numerically find the ground state of 
the system. For antiferromagnetic interactions, opposite 
vortices appear in the m = 1,-1 components with equal 
amplitude when the spin-orbit coupling is weak. Such a 
state carries pure spin current and no mass current. This 
spin current is a super-flow as it is the ground state and 
must be stable. 

We note that there has been a lot of theoretical and ex¬ 
perimental work on the counterflow of two species BEC 
[10H2lj. For two miscible BECs with counterflow, it is 


found that there is a critical relative speed between the 
two species^ beyond which the state is dynamically un¬ 
stable |ll ll9]. It is shown that the instability can lead to 
proliferation of solitons El EH and quantum turbulence 
jl6j |. This kind of counterflow is very similar to a spin 
current but it is not for two reasons: (1) Theoretically, if 
we regard the two species as two components of a pseudo¬ 
spin, this pseudo-spin has no SU(2) rotational symmetry 
as the number of bosons in each species is conserved. 
(2) Experimentally, it is hard to control the number of 
bosons in each component to create a spin current that 
has no mass current. 

The paper is organized as follows. In Sec. [Ill we first 
study the stability of a spin-1 planar counterflow. We 
identify the mechanisms associated with the instabilities, 
and find that the quadratic Zeeman effect can stabilize 
such a planar flow. We then study the similar situation in 
the circular geometry in Sec. Hill The pure spin current 
consists of a vortex and anti-vortex in the m = —1,1 
components. The experimental schemes to realize the 
stable pure spin current is discussed in Sec. IIV1 Finally, 
we briefly summarize our main results in Sec. |Vl 

II. PLANAR FLOW 

The dynamics of a spin-1 BEC in free space is governed 
by the mean field Gross-Pitaevskii (GP) equation [7], 

g fi 2 V 2 \ 

ih—1p m = - 2M~^ m + C°P^ m + C 2 X S'Srnnlpn, (1) 

n =—1 

where {jn = 1,0,—1) are the components of the 
macroscopic wave function, p = ^2 m =-i IV’ml 2 is lh e 
total density, s* = is the spin density 

vector and S = (S x ,S y ,S z ) is the spin operator vector 
with Si (i = x , y , z ) being the three Pauli matrices in 
the spin-1 representation. The collisional interactions in¬ 
clude a spin-independent part Co = 47rft 2 (ao + 2 * 22 ) /3M 
and a spin-dependent part C 2 = 47rfi 2 (a2 — ao)/3M, with 
af(f = 0, 2) being the s -wave scattering length for spin-1 
atoms in the symmetric channel of total spin /. 
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We consider a spin current state of the above GP equa¬ 
tion with the form 


i/j = 




( 2 ) 


where n is the density of the uniform BEC. The require¬ 
ment of equal chemical potential leads to |ki| = |k 2 1 - In 
the case where ki = — k 2 , this state carries a pure spin 
current: the total mass current is zero as it has equal 
mass counterflow while the spin current is nonzero. 

It is instructive to first consider the special case when 
there is no counterflow, i.e., ki = k 2 = 0. The exci¬ 
tation spectra are found to be e° = ^ 2 c 2 ne q + e 2 and 

ef 1 = ^2c§ne q + e q , e^ 1 = ^ 2 c 2 ne q + e^, respectively, 


with e q = h 2 q 2 /2M. So for antiferromagnetic interac¬ 
tion (co > 0 , C 2 > 0 ), all branches of the spectra are real 
and there is a double degeneracy in one branch of the 
spectr a. The phonon excitations give two sound veloci¬ 
ties, y/ncjjM (i = 0 , 2 ), corresponding to the speeds of 
density wave and spin wave, respectively. However, the 
existence of phonon excitation does not mean that the 
pure spin current (ki = k 2 ^ 0 ) is a super-flow as we can 
not obtain the current with ki = k 2 ^ 0 from the state 
with ki = k 2 = 0 by a Galilean transformation. 

The stability of the spin current has been studied in 
Ref. 0 for the case ki = — k 2 ^ 0. It is found that, for 
the antiferromagnetic interaction case (co > 0,C2 > 0 ), 
the excitation spectrum of the m = 0 component always 
has nonzero imaginary part in the long wavelength limit 
as long as there is counterflow between the two compo¬ 
nents, and the imaginary excitations in the m = 1,-1 
components onl y appea r for a large enough relative ve¬ 
locity vi = 2^/nc2/M. For the ferromagnetic interac¬ 
tion case (co > 0 ,C 2 < 0 ), both excitation spectra of the 
m = 0 and m = 1,-1 components have nonzero imagi¬ 
nary parts for any relative velocity. This means that the 
pure spin current cannot be stable in any cases. 

For the general non-collinear case (k = kl + k2 ^ 0) 
and antiferromagnetic interaction, the excitation spec¬ 
trum for the m = 0 component is found to be 


e° = ^ (l k l 2 H k i| 2 )+ c 2™) — c 2 n2 - (3) 

We see here that as long as the momenta of the two com¬ 
ponents are not exactly parallel, i.e., ki is not exactly 
equal to k2, then |k| < |ki|, and there is always dynam¬ 
ical instability for the long wavelength excitations. 

Therefore, the spin current in Eq. m is generally un¬ 
stable and not a super-flow. This instability originates 
from the interaction process described by i n 

the second quantized Hamiltonian. This energetically fa¬ 
vored process converts two particles in the m = 1,-1 
components, respectively, into two stationary particles in 
the m = 0 component. To suppress such a process and 


achieve a stable pure spin current, one can utilize the 
quadratic Zeeman effect. With the quadratic Zeeman ef¬ 
fect of negative coefficient, the Hamiltonian adopts an 
additional term A m 2 (A < 0 and m = 1,0,—1). This 
term does not change the energy of the m = 0 compo¬ 
nent, but lowers the energy of the other two components 
m s 1 ,- 1 . As a result, there arises a barrier for two 
atoms in the m = 1,-1 components scattering to the 
m = 0 component, and the scattering process is thus 
suppressed. 

The above intuitive argument can be made more rigor¬ 
ous and quantitative. Consider the case ki = —k 2 - With 
the quadratic Zeeman term, the excitation spectrum for 
the m = 0 component changes to 


+ C2n — Xj —c 2 n 2 . (4) 

So as long as —A —fi 2 |ki| 2 / 2 M > 0, long wavelength exci¬ 
tations will be stable for the m = 0 component. From the 
excitation of the m = 0 component, one can obtain a crit¬ 
ical relative velocity of the spin current, vq = 2^/—2X/M. 
There is another nonzero critical velocity v\ = 2^/nc2/M 
determined by the excitations of the m = 1,-1 compo¬ 
nents. The overall critical velocity of the system is the 
smaller one of vq and v\. Therefore, below the critical 
relative velocity v c = minj^o, iq}, the pure spin current 
is stable and a super-flow. The experimental scheme to 
realize such a Zeeman effect will be discussed in Sec. m 



III. CIRCULAR FLOW 

In the cylindrical geometry, we consider a pure spin 
current formed by two vortices with opposite circulation 
in the m = 1,-1 components. From similar arguments, 
one can expect that interaction will make such a cur¬ 
rent unstable. Inspired by the quadratic Zeeman effect 
method above, we propose to use spin-orbit coupling to 
stabilize it. The spin-orbit coupling can be viewed as 
a momentum-dependent effective magnetic field that ex¬ 
erts only on the m = 1,-1 components. Therefore, it 
is possible that spin-orbit coupling lowers the energy of 
m = 1,-1 components, and consequently suppresses the 
interaction process leading to the instability. 

The model of spin -1 BEC subject to Rashba spin-orbit 
coupling can be described by the following energy func¬ 
tional, 

T T( ^xPy ~ SyPx) / 5 (5) 


where p is the density, V(r) = \Muj 2 (x 2 + y 2 ) is the trap¬ 
ping potential, and 7 is the strength of spin-orbit cou¬ 
pling. (• • •) is the expectation value taken with respect to 
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the three component wave function ip = (Vh, V’o? V ; -i) T - 
The strength of the spin-orbit coupling 7 defines a char¬ 
acteristic length a soc = h/M^, and can be rescaled to 
be dimensionl ess with respect to the harmonic oscillator 
length ah = yJh/Muj. Then we characterize the strength 
of spin-orbit coupling with the dimensionless quantity 
k = ah/a soc = The spin-orbit coupling of 

Rashba type here can be generated in various ways, which 
will be discussed in the next section. 

The above model can describe a spin -1 BEC of 23 Na 
confined in a two-dimensional harmonic trap. Assume 
the atom number is about 10 6 . Using the estimate of 
scattering lengths ao = 50a#, a 2 = 55a# [22], with a# 
being the Bohr radius, the ground state of spin-1 23 Na 
should be antiferromagnetic because Co > 0 , C 2 > 0 jEj|. 
Previous studies of spin -1 BEC with Rashba spin-orbit 
coupling mostly focus on the strong spin-orbit coupling 
regime, where the ground state is found to be the plane 
wave phase or the stripe phase, for ferromagnetic in¬ 
teraction and antiferromagnetic interaction, respectively 
[23j . Here we are interested in the antiferromagnetic in¬ 
teraction case and the weak spin-orbit coupling regime 
(k « 1 ), and calculate the ground state wave function of 
the energy functional with the method of imaginary time 
evolution. 

We find that when the spin-orbit coupling is weak (/c <C 
1 ), the ground state wave function has the form 

/ xi(r)e"^ \ 

= Xo (r) , ( 6 ) 

\ X-i (r)e^ J 

with Xi( r ) = — X-i(r) and all Xi real. The ground state 
is shown in Fig. [lj Such a ground state consists of an 
anti-vortex in the m = 1 component and a vortex in the 
m = — 1 component. The m = 0 component does not 
carry angular momentum. Since \ipi\ = IV’-il? the net 
mass current vanishes. 

The wave function in Eq. © can be understood in 
the single particle level. In terms of the ladder operators 
of spin and angular momentum, the spin-orbit coupling 
term reads 


'Hsoc = 


7 xMHuj 


S+ (a R - + S- (7 - a L ) 


(7) 


where S± is the ladder operator of spin, and a|^ is the 
creation operator of the left (right) circular quanta [24j . 
When the spin-orbit coupling is very weak (k <C 1), its 
effect can be accounted for in a perturbative way. From 
the ground state = | 0 , 0 ), the first order correction 
to the wave function for small 7 is given by 


' t<,) = 2 ^(- s ^ +s -»«)' o ’ o > 
=|(-ii,-i) + i-i,i)), 


(8) 


where \m s ,m 0 ) denotes a state with spin quantum num¬ 
ber m s and orbital magnetic quantum number m 0 . One 




FIG. 1: (color online) Amplitudes (al,bl,cl) and phase an¬ 
gles (a2,b2,c2) of the three component wave function pj = 
(}pi,'ipo,'ip-i) T at the ground state of Hamiltonian J5} for a 
BEC of 23 Na confined in a 2D harmonic trap. The particle 
number is 10 6 , the frequency of the trap is 2n x 42 Hz, and the 
dimensionless spin-orbit coupling strength is a = 0.04. The 
units of the X and Y axes are ah- 


immediately sees that ^1 has angular momentum —h and 
'ip -1 has angular momentum h. Besides, the amplitudes 
of both 'ipi and ^-1 are proportional to k. 

There exits a continuity equation for spin density and 
spin current, which is 


-(7s^) + v-7 = o. 


(9) 


The spin current density tensor J f, (j i = x,y,z denotes 
the spin component) is defined as |25l . l26j 

K =\ {V’ t £W’ + c.c.} 

= I { E ( S n)mn W nl^l + C.C. > , (10) 


where 


m,n,l 


Vnl = +l(zX S ni) . 


( 11 ) 
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FIG. 2: (color online) Distribution of the spin current densi¬ 
ties JJ (blue arrow), J* (red arrow) and J s z (black arrow) of 
the ground state shown in Fig. |T] The length of the arrows 
represents the strength of the spin current. The arrow length 
of different colors is not to scale, ft = 0.04. The units of the 
X and Y axes are ah- 


and c.c. means the complex conjugate. The second part 
in w n i is induced by the spin-orbit coupling. 

By the definition in Eq. ([TO]) , the spin current density 
carried by the ground state © is 

J* = 7 sin 2 ^>|^i| 2 x + 7 (|V’o| 2 + 2|i/’i | 2 sin 2 </>) y, 

J l = _ 7 (I"001 2 + 2|V’i| 2 cos 2 <p) x — 7sin20|V’i| 2 y, 

J * = (-^T- + V / 27lV J i'0o|) 0- (12) 

From both analytical and numerical results of the wave 
function, |?/q| <C |^o|? so roughly points in the y di¬ 
rection, while J® almost points in the — x direction. J s z 
represents a flow whose amplitude has rotational symme¬ 
try. From the numerical results shown in Fig. [2j we see 
that J z is a counter-clockwise flow. The amplitudes of 
J* and J® are of the same order, both proportional to ft, 
while that of J|, proportional to ft 2 , is much smaller. It 
is evident that the state in Eq. © carries no mass cur¬ 
rent and only pure spin current. Since the spin current 
is in the ground state, it must be stable. In this way, we 
have realized a superfluid of pure spin current, or a pure 
spin super-current. 


IV. EXPERIMENTAL SCHEMES 

In this section, we propose the experimental schemes 
to generate and detect the pure spin currents discussed 
in Sec. Hand Sec. HD 


The planar pure spin current can be easily generated. 
By applying a magnetic field gradient, the two compo¬ 
nents m = 1,-1 will be accelerated in opposite direc¬ 
tions and _a pure spin current is generated as done in 
Refs. [14, 15]. To stabilize this spin current, one needs to 
generate the quadratic Zeeman effect. We apply an oscil¬ 
lating magnetic field B sin ujt with the frequency uj being 
much larger than the characteristic frequency of the con¬ 
densate, e.g., the chemical potential y. The time averag¬ 
ing removes the linear Zeeman effect; only the quadratic 
Zeeman effect remains. The coefficient of the quadratic 
Zeeman effect from the second-order perturbation theory 
is given by A = {gy^B) 2 /AFhf, where g is the Lande g- 
factor of the atom, /ip is the Bohr magneton, and AE^f is 
the hyperfine energy splitting [27]. For the F = 2 mani¬ 
fold of 87 Rb, AFhf < 0, so the coefficient of the quadratic 
Zeeman effect is negative. 

The circular flow in Sec. mi may find prospective real¬ 
izations in two different systems: cold atoms and exciton 
BEC. In cold atoms, we consider a system consisting of a 
BEC of 23 Na confined in a pancake trap, where the con¬ 
finement in the z direction is so tight that one can treat 
the system effectively as two dimensional. The spin-orbit 
coupling can be induced by two different methods. One is 
by the exertion of a strong external electric field E in the 
z direction. Due to the relativistic effect, the magnetic 
moment of the atom will experience a weak spin-orbit 
coupling, where the strength 7 = #/tb|E| /Me 2 . Here M 
is the atomic mass and c is the speed of light. For weak 
spin-orbit coupling (small 7 ), the fraction of atoms in the 
m — 1,-1 components is proportional to 7 2 . For an ex¬ 
perimentally observable fraction of atoms, e.g., 0 . 1 % of 
10 6 atoms, using the typical parameters of 23 Na BEC, 
the estimated electric field is of the same order of mag¬ 
nitude as the vacuum breakdown field. For atoms with 
smaller mass or larger magnetic moment, the required 
electric field can be lowered. Another method of realizing 
spin-orbit coupling is to exploit the atom laser interac¬ 
tion, where strong spin-orbit coupling can be created in 
principle [Hj. In exciton BEC systems, as the effective 
mass of exciton is much smaller than that of atom, the 
required electric field is four to five orders of magnitude 
smaller, which is quite feasible in experiments |29 h32| . 

The vortex and anti-vortex in the m = 1,-1 compo¬ 
nents can be detected by the method of time of flight. 
First one can split the three spin components with the 
Stern-Gerlach effect. The appearance of vortex or anti¬ 
vortex in the m = 1,-1 components is signaled by a 
ring structure in the time of flight image. After a suffi¬ 
ciently long time of expansion, the ring structure should 
be clearly visible. 


V. CONCLUSION 

In summary, we have studied the stability of a pure 
spin current of a spin -1 BEC. In the planar flow, the 
system always suffers from dynamical instability. The 
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origin of the instability is the interaction process that 
converts two particles in the m = 1,-1 components into 
the m = 0 component. Based on this, we propose a 
method to stabilize the pure spin current by utilizing the 
quadratic Zeeman effect. In the circular flow, we have 
proposed to use spin-orbit coupling to make the pure spin 
current stable. For weak spin-orbit coupling, we have 
found that the ground state of the system is a superfluid 
of pure spin current. The experimental schemes to realize 


and detect these pure spin currents have been discussed. 
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